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Abstract. The asymptotic behavior of the micromagnetic free energy gov- 
erning a ferromagnetic film is studied as its thickness gets smaller and smaller 
compared to its cross section. Here the static Ma^cwell equations are treated as 
a Murat's constant-rank PDE constraint on the energy functional. In contrast 
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field without solving the magnetostatic equations. In particular, the math- 
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213-223] regarding convergence of minimizcrs are recovered by giving a char- 
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1. Introduction 

Over the last twenty years there has been tremendous scientific progress in the 
research on thin-film devices pushing forward technology and leading to important 
industrial applications. By a thin film one understands a layer of material whose 
thickness ranges between a fractional amount of a nanometer and a couple of mi- 
crometers. When using thin films in computer data storage media and solar cells, a 
deep understanding of their ferromagnetic properties is of great importance [151124] . 
A widely used mathematical procedure to achieve exactly that takes the theory of 
micromagnetics for bulk bodies [3l|4] as a starting point and derives reduced the- 
ories capable of capturing the specific features of thin material layers by means of 
dimension reduction techniques. 

It was in |13[ that for the first time authors studied convergence of minimizers 
of the micromagnetic energy on a film whose thickness gets smaller and smaller 
relative to its cross section. In this article we follow the same approach regarding 
scaling, which corresponds to considering the regime of very small film samples, but 
we employ an equivalent formulation in the sense of [S]. The latter illustrates that 
micromagnetism is actually one of the examples, where the mathematical modeling 
of physical phenomena within a variational formulation leads to functionals that 
do not simply depend on gradient fields. Instead one faces more intricate partial 
differential constraints that involve an interaction of divergence- and curl- free vector 
fields. Precisely, the static Maxwell (or magnetostatic) equations, which govern the 
relation between the magnetization m of a ferromagnetic body occupying a bounded 
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domain C M"^ and its induced magnetic field /i : M'' — ^ M"^, read 

div(TO + ft.)=0 inM^ (1.1) 

curl ft. = in]R^ 

where m : £7 — > M'^ is identified with its trivial extension to the whole space by zero. 
In the literature this mathematical difficulty is commonly tackled by explicitly 
solving or proving existence of solutions to the magnetostatic equations in their 
weak formulation and expressing h in terms of to. Our idea is to keep both to 
and h as variables and work directly with the PDE constraint imposed by (|l.ll) . 
The main result of this paper is a rigorous F-convergence based 3d-2d dimension 
reduction for the constrained micromagnetic functional. Our approach allows us to 
keep track during the limit process not only of the magnetization, but at the same 
time also of the induced field. 



2. Formulation of the problem and statement of the main result 

Let fig = a; X (0,£) model a ferromagnetic body of thickness e > with cross 
section w C M^, wlog \uj\ = 1. The free energy per unit volume that emerges in the 
theory of micromagnetism [3J[2S] is given by 

p r- 71 /-/ a\\7m\^ + ipifh)dy+l [ \h\^ dy, if (m.h) & Ve, 
[to, h\ = i e Jfj^ 2 7r3 

I oo, otherwise. 

Here a > is a material constant and (y9 : M'^ — > [0,cx)) is a continuous, even 
function featuring crystallographic symmetry. Further, 



|(m,ft) e X L2(]r3.j^3) 



_^mag/ rn ) ^ in KMtoI = ?7i, in a 



h 

The first order PDE constraint in the definition of V^, namely 

fh \ _ f div div \ / rh \ _ / div(TO + h) 
h I l^O cmT) [ h I ^ [ curl ft 



A" 



0, 



conveys the magnetostatic equations. The operator curl is supposed to be inter- 
preted as curl = Vx , i.e. 

curl ft = (i92ft3 - d3h2, d^hi - dihs, dih2 ~ (?2fti)'^. 

Notice that throughout this work all occurring differential operators and partial 
derivatives are to be understood in the sense of distributions, for example div(m + 
ft) = in means - J^-, (to + ft) • V(/i dy = for aU test functions (j) G C^i^^). 

Physically speaking, the nonconvex constraint |to| = rus in 0,^ (depending on the 
regularity of to, this equality may only be fulfilled pointwise a.e. in f2g) encodes the 
fundamental assumption that the body is locally saturated with saturation magne- 
tization TOs > 0. The second term in the definition of is the anisotropy energy, 
which penalizes magnetizations varying from special directions within the crystal 
lattice of the ferromagnet. The latter are called the easy axes of magnetization. 
The contribution of exchange energy is captured by the first term of Eg. It results 
from a force tending to align magnetic moments of neighboring atoms and therefore 
favors regions of constant magnetization. The third summand in Eg is an integral 
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over the whole space M.^ modehng the energy of the magnetic field h induced by 
TO. These three energy components impose competing requirements on the magne- 
tization and minimizers of may form interesting structures. In the case of bulk 
bodies one observes Weiss domains separated by thin Bloch walls. For more details 
on the physical motivation and interpretation of the nonlocal and nonconvex energy 
Eg see for example [Tl[T6l[l7j[26] and the references therein. 

As originally stated in [25] and further discussed in [12], jg ^ first order 

differential operator that meets Murat's constant-rank property [55], i.e. the symbol 
A^mag of y^^^^g satisfies 

rank A^mag (^) = const. for aU ^ £ (2.1) 
Indeed, A^mag(^) e Lm{R^ x M x E^) and it holds 

ker A^mag (^) = { (a, &) e X : ^ • (a + &) = 0, ^ x 6 = } 

= { (a, 6) e X . 5 = AC, C • a = -A, A e E } . 

So rankA^mag(C) = 6 — dimker A_4mag(C) = 3 for all ^ S S^. Hence, the problem 
we are interested in can be studied in the more abstract context of dimension 
reduction for functionals on ^-free vector fields. The work on variational problems 
within the ^-free framework can be traced back to [S] and was advanced by Fonseca 
and Miiller |12j . who came up with the notion of ^-quasiconvexity (in its modern 
sense) and studied lower semicontinuity of functionals involving integrands with 
this property. Since then, a lot of papers investigating for instance relaxation, 
homogenization and Young measures in the y^-free setting have been published [2l 
[T0HT2] . The first article to cover 3d-2d asymptotic analysis in such generality 
is |18j , while the special case of thin- film limits for gradient dependent problems 
{A = curl) has been treated before (see for instance [HlllIllHI). For a recent result 
in the context of functionals on solenoidal vector fields {A = div) we refer to [TO] . 
In fact, [18| provides the technical basis for the work presented in the following. 

To obtain a variational problem on the fixed domain fii = uj x (0,1), we apply 
the standard parameter rescaling, 

X = {x',Xd) = {y\e~^yd) 

for y = (yi,^ .. ,2/^-1, yrf) = {v' ,Vd) e Vt^, and we set m{x) = m{y) = m{x' ,exd) 
and h{x) = h{y) = h(x' , sxd)- Then transforms into 

J / a iVeTOp + </7(to) da: + i / \h\'^ dx, if (to, /i) e Z^e, 
[ 00, otherwise , 

where 

^^^s ^ = in |77i| TO, in f^i}. 

The rescaled versions of the operators V = (V, 83)'^ = {di, 82, 83)'^ and A™'^^ are 
given by 

Ve = (V, i^a)^ and A^^"^ = A^mag(Ve), (2.2) 
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respectively. Accordingly, we define 

dive = Ve • and cui% = x . 
In view of these definitions the main result is the following: 

Theorem 2.1. The T-limit of for e — > 0^ with respect to weak convergence in 
W^''^{fli;R'^) X L^(R-^;R^) exists and is represented as 



Fo[m,h] = 



I a |V'mp + (/j(to) da; + i ( \h\^ dx, if{m,h)eUo, 



otherwise , 



with 



A^^^ ^ ^ = m M'^, d^m = in fti, \m\ ~ nig in f2i| 



and J^^^ defined through 



4" 

•^0 



mag 



/ 93(^3 + /^s) \ 

-^3/12 

V dih2 - d2hi J 



Moreover, compactness holds in the weak topology of W^''^{ili;M.^) x L-^(M 

As we discuss in Section [5] the result of Theorem 12.11 is in complete agreement 
with the limiting micromagnetic energy derived in |13j . 

In [18j . dimension reduction within the general y^-free setting is investigated. 
One purpose of Theorem 12.11 is to illustrate the power of the concepts and tools 
developed there by means of another physically relevant example. (An immediate 
first application to bending of thin films in nonlinear elasticity is studied in |18[ 
Section 5].) 

Since the functionals contain first order derivatives of m and involve a non- 
convex constraint, they do not exactly fit into the context of [T5]. (The issue of 
i^e being defined on functions on the whole space can be overcome by replacing 
Fourier series with Fourier transforms in the proofs. However, it is not clear how to 
extend that theory to mixed-order differential operators.) Therefore, we consider 
i^e as split into a part that is convex in the derivatives of m and one that is in 
line with [TS]. When it comes to proving the upper bound, the crucial step is to 
exploit a tool introduced in [18], which yields convergence of the symbols of A^'^^ 
for e — > 0"^. This provides a "first candidate" for the recovery sequence which, 
however, lacks the necessary regularity and fails to meet the nonconvex constraint 
imposed by the requirement of local saturation. To handle this matter, we modify 
the sequence by choosing the magnetizations to be constant and by adjusting the 
exterior fields with the help of projection operators onto curlg-free fields. 

In the next two sections we give the detailed proof of Theorem 12.11 by showing 
separately the required upper and lower bounds. We remark that throughout this 
work we use generalized sequences with index e > 0, like (ue)e, by which we refer 
to any sequence (we^Oj with ej — > 0+ as j — > 00. 
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3. Proof of compactness and the lower bound 

We begin by proving the following compactness result, which is essentially based 
on the coercivity of the micromagnetic free energy. Notice that extracted subse- 
quences are not relabeled in the sequel. 

Proposition 3.1 (Compactness). Let Sj — > 0+ for j — > oo. Further assume 
that {m^.,h^.)j C W^''^{fli;W^) x i^(M''';M^) is a bounded energy sequence for F^. , 
i.e. 

Pej [^Ej ,h^^] < C < oo for all j G N. 

Then there exists a subsequence {m^. , h^.)j and (m, h) G iy^'^(rii; R"^) x L^(]R^; R"^) 
such that 

m^^^m inW^''^{Q.i;W^), (3.1) 

for J — >■ oo. Moreover, it holds that (to, h) (z Uo- 

Proof. In view of the constraint \m^j \ = in fli and the fact that 

||Ve^TOeJjL2(ai;B3x3) < C < c» foralljeN, (3.2) 

one infers that Uto^^ ||vi/i,2(nj;R3x3) is bounded uniformly with respect to j, which 
implies the existence of a subsequence {m^.)j and a function to S W^^'^(fii; R^) 
such that m^j ^ m in VF^'^(r2i; R^). By compact embedding we find (after passing 
to a subsequence) that to^^ — >■ m pointwise a.e. in ili as j — oo, so that |m| = to^ in 
ill. Recalling the definition of Ve in (|2.2p . we conclude from p.2p that dsrue^ — > 
in i^(J7i;R^). Thus, dsm = in f2i by uniqueness of the limit. Since the induced 
energy contribution of F^^ [m^. , h^. ] is bounded, one can extract a subsequence of 
{hej)j satisfying 

for some h G L^(51i; M?). The expression curl^^ h^. = in R'^ is equivalent to 

92(/le, )3-l/£j53(/le, )2 =0 in R^, 

l/e,a3(/ie,)i-ai(/i,j3 = inR3, 
di{K^)2-d2{K^)i^Q inR3. 

When passing to the limit j — > oo, it follows that 53/12 = c^a/ii = and dih2 = d2hi 
in R'^. Finally we exploit div^^ (m^^. + /le^ ) = in R'^ to derive 

53(TO3 + /i3)=0 inR^ 

Thus, (m,/i) G Z^o- □ 

Remark 3.2. If (we^. , )ej is a sequence of minimizers for F^^ , the convergence 
in (|3.ip can be shown to be strong in W^''^ and L^, respectively (compare [T51 
Theorem 4.1]). 
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Next we give the proof of the lower bound. Let nie — ^ m in VF^'^(f2i; M'^) and 
he h in L^(R'^;R'^) as e tends to zero. Consider any ej — > 0+ for j — > oo and 
assume that 

Hminf i^e [TOe /le ] = hm i^^ [toj /i^J < oo, 

otherwise the corresponding hminf-inequaUty is immediate. Then (jUe^ , he^ )j C 
is of bounded energy and (to, h) € Uq by Proposition 13.11 In view of the compact 
embedding W^''^{V.i;M.^) ^ L^(f2i;R'^) we obtain a subsequence with rriej to 
pointwise a.e. in fJi, so that by the continuity of ip, the weak lower semicontinuity 
of the L^-norm and Fatou's lemma, 



lim jTOe , /le J > lim / a IV'to^ ^ + <y5(TOe J dx + i / {h^A'^ dx 
>/ a iV'TOp + (^(m) da; + i / da;. 



Hence, 



liminf i<"^[TO^, hJ > FoWn, h], 



which is the liminf-inequality. 

4. Construction of a recovery sequence 

This section is based on arguments involving operators of the following form: For 
given matrices . . . , A^'^) G M'x" let A be the linear partial differential operator 
of first order defined through 

d 

^:=^AW9fe. (4.1) 

fc=i 

Then the symbol A_4 of A is given by 

d 
k=l 

The essential assumption on A is Murat's constant-rank condition P^ll^ . precisely 

rank A aH) = const. for aU ^ G S'^"^ 

By [T51 Lemma 2.2] the operators Ae '■= A^(Ve) are of constant rank for all e > 
provided A has the same property. 

As established in Section [2] (see (|2.ip ) A"^'^^ fits into the framework described 
above with 



d = 3, 71 = 6, 1 = 4: and u = 



TO 

h 



So, after having proved some technical tools for general constant-rank operators 
we will be able to apply these findings to the context of micromagnetics. 

The next theorem is a modification of [HI Theorem 2.7], which is formulated 
for LP-functions on the torus, for L^-functions on the whole space. The important 
issue in comparison with |12l Lemma 2.14] is to obtain constants independent of e. 
A first comment in this direction is made in [19j . where projection operators onto 
div-free fields are investigated. 
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Lemma 4.1 (Projection onto ^e-free fields). Suppose A is a constant-rank 
operator as defined in (j4.ip and s € (0,1). Then there exist bounded operators 
Va, ■■ L2(R'^;R") L2(M'^;R") with the following properties: 
(i) {Va. ° VaJu = Va^u for all u e L2(]Rd;M«). 
(a) [Ae o Va,)u = Q inW^ for all u £ L2(R'^;R"). 
(Hi) The operators Vac uniformly bounded with respect to e, i.e. 

for all u e L^^Rd.jgri) ^j^^^ 

a constant C > independent of e. 
(iv) There exists a constant C > such that 

11"" ~ '^'4E"IL2(Rd;R") - ^ ll-^e""!! W-1.2(Rd;Ri) 

for all u <E L2(R'^;R") and all e S (0, 1). 

Proof. In the sequel we employ the common notation J- to refer to the Fourier 
transform and use J^~^ for its inversion. 

For €W^\ {0} let Pxe(C) be the orthogonal projector onto kcr A^^(^) C R". 
So the mapping P^^ ; M'' \ {0} Lin(R";R") is 0-homogeneous and smooth. The 
fact that the operator norm of P^^ (^) is equal to 1 for all ^ S R*^ \ {0} and e > 
renders Va, ■ i^(K'';M") L'^{R'^;W) defined by 

7'^.u = ^-i(P^,(.)-Fu), uei2(R'';R"), (4.2) 

a continuous operator satisfying the estimate 

for all u G L^(R'';K"). This proves (Hi). The properties (i) and (ii) are an 
immediate consequence of the structure of P^^ together with (|4.2I) . 

In order to show {iv) let QaAO & Lin(R';K") with ^ e R'' \ {0} and e > be 
given by 



^aAO^ for V G range A.A^(^) with v = Aa^I^z, z € 
for V e (range A.Ae(0) • 



Then, Qa, : R'' \ {0} Lin(R';M") is homo geneous of degree —1 and smooth. 
Notice that the smoothness of both Qa, and Pa, rests upon the constant-rank 
property of Ae (compare [HI Proposition 2.7]). Besides, Q^^(VM) is bounded 
in the L°°-norm uniformly with respect to e. This can be seen as follows. Since 
^aAv) = ^AiVe) for aU ?7 e R'' with 77^ := (771, 77^-1, l/eTy^)^, we may argue 
for e e (0, 1) that 

sup ||Q.A,('7)|lLin(R';R") = SUp |lQ.A('7£)||Lin(R';R") 

< sup ||Q^(C)llLin(R';R") = sup a"lQ.A(?/)||Lin(R';R") 

4eR'^,l€l>i r)eS''-i,Q>i 

< sup ||Q^(?7)||Li„(Ri:R..) < 00, 



8 



CAROLIN KREISBECK 



where the second equaUty is a consequence of the (— l)-homogeneity of Qa- The 
final estimate resuhs from Qa. being smooth on the unit sphere which is a 

compact subset of M.'^. 

Consequently, for <E L^(]R''; R') given through J^w^ = |'|~"^A^^ it holds 

that 



The last inequality holds by the definition of the W "'^'^-norm. Using the properties 
of Qas ; ^Ae ^iid the linearity of A_4^ yields 



Remark 4.2. For p G (1,cxd) the statement of Lemma [4. II is still true. However, 
the line of reasoning is more involved due to the fact that the Fourier inversion 
formula docs not hold in general. In analogy to the proof of p[8^ Theorem 2.7] 
(with Fourier series replaced by Fourier transforms, where necessary), this difficulty 
can be overcome by using approximation via smooth functions in conjunction with 
Mihlin's Multiplier Theorem (see for instance [111 Theorem 5.2.7]) and a scaling 
argument for Fourier multipliers. 

As an essential tool towards the construction of a recovery sequence, we show 
the following analog of [HI Proposition 4.1] within the setting of functions defined 
on the whole space. Naturally, an extension property in the sense of [181 Assump- 
tion A3] is not needed. The proof follows closely along the lines of [18], but is 
substantially easier, since one of the two relevant terms is forced to vanish here. 

Proposition 4.3. Let A be a constant-rank operator as in (|4.ip such that the 
number of non-zero rows of the matrix A^'^' is equal to the rank of A^'^'^ . Further 
suppose u e L^(R'';R") satisfies Aou = in R'', where 




(4.3) 



< C\\AeU 



= J^u-PaA-ZI-D-^u^ J'u-PaA-)-^'^- 
In view of (|4.3p and (|4.2p this proves (iv). 



□ 




and [A^'^^Y denotes the ith row of A^'^K Then there exists a sequence {u^)^ C 
L2(R'';R") with AeUe = m and ^ u m L'^{W^:W) for e 0+. 



For the proof of this proposition the auxiliary "symbol 
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will be needed. By Po(^) we denote the orthogonal projection onto ker Ao(^) C R". 
The symbol A_Ao coincides with Aq outside the hyperplane where = 0, in formulas 

AoiO = (0 for aU ^ e R'^ with ^ 0. (4.4) 

Note that in contrast to A^^, the expression Aq is not the symbol of a constant- 
coefficient partial differential operator (see [THl Remark 4.2]). 

Actually, Aq turns out to characterize the limit behavior of the symbols A^^ as e 
tends to zero. The exact result is formulated in the next lemma, which was proven 
in [18] and is repeated here for the readers' convenience. 



Lemma 4.4 ( [181, Lemma 4.3] ). If A meets the assumptions of Proposition^ 
the symbols A^^ converge to Aq as e — > 0+ in the sense that ^A^iO^ ^ ^o{S,)^' 
all ^ e R'' and all v G R". 



Proof of Proposition \4-3\ The fact that u is ^o-frec implies 

Aao{0{J'u){0 ^0 for a.e. e e R''. (4.5) 

Now we split Tu into 

{Tuf\0:-i^^''^^^^ ''^'^'^ and {Tuf\0--i' ''^'^'^ 
^ ' [O ife<i = 0, ^ ^ 1(-F«)(e) if^d^O, 

so that by means of Fourier inversion, 

With this definition, u^^^ and w^^) are L^^jjd- ]^™)_f^j^(2|;jQjjs satisfying J"w(*) — 
(J"w)''^ for i = 1,2, and we may conclude from ()4.5p that 

^0^(1) = ^0^(2) =0 in R''. 

Turning to u(2) we observe that by construction ddU^'^^ = 0. In order to have 
the quadratic integrability of u^^) on R'' preserved it needs to hold that u'-^) = 0. 
Hence, for the proof of this proposition it will be enough to show the existence of 
an ^e-free sequence (ue)e with — > u^^^ in L2(R'*;R") as e — > 0+. 

For this purpose we set := T^.a^m'^' for every e > or, speaking in terms of 
Fourier transforms, 

Fu, :-P^^(.)-Fm(i) 

with and P^^ as in Lemma [4.11 In view of (|4.4p and ()4.5p one finds for a.e. 
S^eW^ that Ao(0(-7^u(^')(0 = 0> which implies 

Recalling that F (and F^^) are L2_igometries we may argue that 

|(P^.(0-Po(C))(-^u(i))(Orcie (4.6) 

At this point we apply Lemma l4.4l This, together with the uniform boundedness of 
the projection operators P^^, allows us to use Lcbesgue's Dominated Convergence 
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Theorem and we conclude that the right-hand side in (|4.6p tends to zero for e — > 
0+. □ 

As mentioned at the beginning of this section, A^^^ is of the form (|4.ip meeting 
the constant-rank condition and one can check easily that both the rank of (A'^'^s^Cs) 
and its number of non-zero rows is three. Besides, a straightforward calculation 
shows that ^q"''^ as in the statement of Theorem 12.11 corresponds to A™^^ in the 
sense of Proposition 14.31 Hence, the next corollary is an immediate consequence of 
the previous proposition. 

Corollary 4.5. For every (m, h) € Uq there exists a sequence 

such that {m^,h^) {m,h) in L'^{R^;R^) x ^^(R^.jgS) as e ^ 0+ and 

for all e > 0. 

Before we prove the upper bound, we establish the following relation between 
curlg- and divg-free fields, which is a type of Hclmholtz decomposition for L^- 
functions defined on M^. 



Lemma 4.6. For every e > it holds that 

where T'curie , T'dive ■ -^^(K"^; K'^) — > i^(R'^;R'^) are projection operators as defined in 
Lemma l4Tl\ and I is the identity map on L^(M^;M'^). 

Proof. In view of (|4.2p . the claim holds true if 

I - Pcuru (0 - Pdiv. (0 for all C e \ {0}. 

with I : R'^ R"^ the identity function. In the sequel we use the notation ~ 
(6,6, for C e Let ^ e R^ \ {0} be fixed. Then 

kcrAdiv.(e) = {i'GK' : 6 • « = } = span{e,^-i, ^^^-^j 
with orthogonal unit vectors ^^'^ and £,^''^, while 

ker Acuru (6 = { « e R3 : ^,xv^O}= span{ee/l6|}- 
Thus, for w e one may infer 

(I - Vcu.lAO)v = V-{V 6)6/|6P = {V ■ i^-'K^ + {V ■ 6^'2)6^'^ = Pdiv.(0^ 

and the proof is complete. □ 

Let us point out that the functions rhe in Corollary 14.51 do not have W^''^- 
regularity nor do they fulfill the required nonconvex constraint of local saturation. 
Consequently, (toe, hg)/; fails to be a correct recovery sequence. In order to overcome 
this problem we make a construction based on the use of appropriate projection 
operators and prove the following proposition. 
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Proposition 4.7 (Recovery sequence). For every {m,h) E Uq there exists 
a sequence {me,he)e C Ue with [me, he) (m,h) in W'^''^{p.i;M?) x L2(K^;M3) 
satisfying 

lim _F£[m£,/i£] = Fi:)[m,h]. 

e->0+ 

Proof. For given (?ti, /i) e Z//o let (to^, /le)^ be as in Corollary 14.51 We set for e > 0, 

= TO, 

he = Vcurl, {he - TO + TOg) . 

The assertion is that this definition of {me,he)e provides a recovery sequence for 
(to, h). 

Indeed, it holds that A'^^^{me, he) = in for every e > 0, since curl^ he ^ 
by Lemma l4.1l fii) and 

diVe(TO£ + he) = diVe (to + PcuiL {he — + ™e)) 

= diVe(TOj + /le) + diVj((7'curl, - I){he - m + me)) = 

in view of Lemma 14.61 and Corollary 14.51 Trivially, satisfies Ito^I = |to| = TOs in 
r^i for all e > 0. Moreover, applying Lemma [4. 1 l (iii)- fiv) yields 

Ij^e - ^!Il2(M3;R3) < \\he — /le||i2(R3.R3) + \\he - L2(r3.R3) 

< ll'PcurU^e — /ie||L2(R3.R3) + \\Vcui%{me - JTl) || i2(R3.R3) + \\he — /l||i2(K3.R3) 

< C (II curie ^e||H'-1.2(K3;R3) + ||?Tie - "l||i2(R3.R3)) + \\he - /l|| l2(r3.r3) 
= C Wrhe - TO||i2(I{3.R3) + ll^-e - /i||l2(r3.r3). 

This expression tends to zero by Corollary 14.51 as e 0+. 
Summarizing, we find that {m.e, he) £ Ue for all e > and 



lim / a iVeTOeP + Gsfme) dx + i / \he\'^ dx 



a |V'toP + (^(to) dx + i / |/ip dec. 



□ 



5. Comparison with the results by Gioia and James 

To demonstrate the agreement between Theorem 12.11 and the result obtained 
in [13l Theorem 4.1] we will prove that the F-limit of i^e as e — > 0+ can be expressed 
equivalently in the form 

1 

—1 
2 

00, otherwise. 



Fo[mh] = { J " l^'^P + + dx, if (to, /i) G ZYo, 



where 

iio ^ {{m,h) eW^'^{uj;M.^) X L^{R^;R^) : |to| = in w, /i = -(0, 0, TOg)^ } 



12 



CAROLIN KREISBECK 



Indeed, if we identify m and h with their constant extensions in space direction X3, 
it is sufficient to show that 

This cquaUty results essentially from the observation that the only function that 
is constant with respect to one of the coordinate directions and at the same time 
quadratically intcgrable on the whole space is the zero mapping. So for (m, h) G Uq 
one finds 93/12 = ft/ii = in M.^ and d^h^ = in M'^ \ ili, which results in 
hi = h2 = and /13 = in R"^ \ fii. respectively. Applying the same argument once 
again to ^((ma + h3)xQi) = in M.^, where is the indicator function of {li, 
yields ~ ^n^sXni = ^'713. Finally, since m is independent of the xa-variable, 
(to, h) is actually two-dimensional and we conclude (to, h) £ Uq. The reverse inclu- 
sion follows simply from the above mentioned identification and a straightforward 
calculation. 

Remark 5.1. a) In view of its representation Fq, the F-limit of is purely two- 
dimensional. The third additive term in the density of _fo indicates that magnetiza- 
tions pointing out of plane are penalized and hence less favorable when it comes to 
energy minimization. Interestingly, we also observe that the limit functional _Fo can 
be viewed as local. In fact, the Maxwell equations disappear in the limit e — > 0+, so 
that the asymptotic problem is free of any magnetostatic constraints of this form. 
Instead, the relation between to and h is ruled by the simple pointwisc equality 
h = — (0,0,TO3)'^. For a detailed interpretation of our F-convergence result with 
respect to physical and engineering applications we refer to |13[ Section 5]. There, 
Gioia and James analyze the qualitative effect of external fields and the practical 
scope of the theory regarding the thickness of films by giving precise estimates for 
some relevant materials. 

b) The main advantage of our approach is the convenient access to complete 
information about the induced field h, which is automatically included in Fq. In 
contrast, the reasoning in [13] requires explicit solving of the magnetostatic equa- 
tions together with a limit analysis for the solutions as tends to zero to gain the 
same insight. 

c) The scaling of the micromagnctic free energy in this work rests fundamentally 
on the assumption that the material parameter a does not depend on the thickness 
e. We refer to |8] for a review paper discussing a number of different reduction 
regimes. Here the characteristic length scale d (~ -y/a) of th*^ magnetic material 
is related to e and I (the length scale of the cross section) in such a way that 
(P/le — )- and e/l — >■ 0. In [7], for instance, the authors study the limit behavior 
in the regime le/(P 3> ln(Z/e) and are able to capture domain structures observed 
in large thin films, as comparison with experimental data underlines. Let us point 
out once again that the results obtained in [13j and recovered here are physically 
relevant only for thin-film samples of sufficiently small lateral expansion. In this 
case one observes a single region of uniform magnetization. 
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